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•  Plastic  shakedown  reduces  stress  and  enhances  capacity  with  good  safety. 
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•  Design  schemes  are  provided  for  plastic  shakedown  and  cyclic  plasticity. 
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This  paper  proposes  a  strategy  to  reduce  the  diffusion  induced  stress  and  enhance  the  capacity  of  a 
layered  electrode  by  allowing  the  plastic  deformation  of  current  collector.  Based  on  analytical  formu¬ 
lations  of  the  stress  in  whole  electrode,  three  types  of  elastoplastic  behaviors  of  current  collector,  i.e.  pure 
elastic  deformation,  plastic  shakedown  and  cyclic  plasticity,  are  identified.  Criterions  separating  the  three 
cases  are  proposed.  It  is  found  applying  a  thin  current  collector  and  allowing  it  to  plastically  yield  in  the 
charge/discharge  cycles  is  beneficial  not  only  to  capacity  as  more  space  can  be  provided  for  active  ma¬ 
terials  but  also  to  electrochemical  stability  because  the  stress  in  active  layer  is  significantly  reduced. 
Structural  design  corresponding  to  plastic  shakedown  shows  good  balance  between  the  said  improve¬ 
ments  and  structural  safety,  whereas  the  case  of  cyclic  plasticity  further  enhances  the  improvements. 
Therefore,  structural  designing  scheme  is  provided  for  the  former  case  according  to  the  criterion  of 
plastic  shakedown  but  for  the  latter  one  based  on  the  Coffin— Manson  relation  with  expected  cycle  life. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  batteries  have  been  widely  used  in  the  applications 
from  portable  electronics  to  electric  vehicles.  The  electrode  of 
lithium  ion  battery  usually  employs  a  multilayer  structure  which  is 
composed  of  active  materials  such  as  graphite,  silicon  or  LiCoC^, 
and  current  collector  made  of  copper  or  aluminum.  In  the  cycles  of 
charge  and  discharge,  stress  field  is  developed  in  the  whole  elec¬ 
trode  partially  because  the  lithiation  induced  deformation  of  active 
layer  is  restricted  by  the  current  collector,  leading  to  mechanical 
and  further  electrochemical  degradation  of  the  electrodes  [1]. 
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Among  the  intensive  studies  devoted  to  the  diffusion  induced 
stress  (DIS),  the  current  collector  was  customarily  considered  as  an 
elastic  material  [2-5].  Mechanical  plastic  yield  is  usually  avoided 
perhaps  because  of  the  traditional  believe  in  mechanical  design 
that  plastic  deformation  is  harmful  to  structural  safety.  However, 
the  main  function  of  current  collector  in  batteries  is  to  provide  a 
passageway  for  electron  transportation,  not  for  load  bearing.  In 
addition,  if  the  current  collector  plastically  yields  in  charge  and 
discharge,  the  relaxation  brought  by  the  plastic  deformation  may 
reduce  the  stress  in  active  layer  from  which  the  electrochemical 
performance  may  benefit.  Therefore,  this  article  attempts  to 
explore  the  possibility  of  allowing  the  plastic  deformation  of  cur¬ 
rent  collector  in  designing  and  taking  advantage  of  the  elastoplastic 
behavior  of  current  collector  to  improve  the  performance  of  a 
battery.  In  literature,  although  the  plastic  yield  of  active  material 
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has  been  widely  discussed  [6-14],  the  impacts  of  plastic  yield  of 
current  collector  have  been  rarely  investigated. 

In  this  manuscript,  we  will  demonstrate  a  strategy  to  reduce  the 
DIS  and  enhance  the  capacity  of  a  layered  electrode  by  applying  a 
thin  current  collector  whose  plastic  deformation  is  utilized.  There 
are  three  types  of  possible  elastoplastic  behaviors  for  the  current 
collector,  i.e.  pure  elastic  deformation  (PED),  plastic  shakedown 
(PS)  and  cyclic  plasticity  (CP).  See  Fig.  1,  taking  a  cyclic  uniaxial 
tension/compression  test  as  example,  pure  elastic  deformation  is 
one  in  which  plastic  yield  does  not  occur.  Plastic  shakedown  is  one 
in  which  plastic  yield  takes  place  only  during  the  first  loading, 
while  the  subsequent  unloading  and  reloading  are  perfectly  elastic. 
Cyclic  plasticity  is  one  in  which  alternating  plastic  yield  of  the 
material  takes  place  in  each  loading-unloading  cycle,  leading  to  a 
closed  elastic— plastic  loop. 

In  order  to  characterize  the  impacts  of  the  elastoplastic  behavior 
of  current  collector  on  battery  performance  and  designing  insights, 
an  analytical  model  will  be  established  for  a  planar  layered  elec¬ 
trode.  Yield  criterion  of  current  collector  will  be  proposed,  taking 
into  account  the  impacts  of  material  properties,  electrode  config¬ 
uration  and  charging  states.  Evolutions  of  the  stresses  in  both 
current  collector  and  active  layer  in  the  cycles  of  charge  and 
discharge  will  be  simulated.  Finally,  designing  insights  will  be 
provided  based  on  the  discussions. 


2.  Methodology 


Consider  a  layered  electrode  in  which  two  active  layers  of  equal 
thickness  h\  are  symmetrically  bonded  to  a  current  collector  of 
thickness  hc,  Fig.  2.  Let  the  thickness  direction  be  aligned  with  the  z- 
axis  and  the  in-plane  of  plate  with  x-  and  y-axes.  Lithium  ions  are 
inserted  into  and  extracted  out  of  the  electrode  from  both  side 
surfaces. 

The  Li-ion  diffusion  in  the  active  layer  is  assumed  to  be  governed 
by  Fields  law.  Other  factors  affecting  the  diffusion  such  as  saturation 
cap  [15],  stress  coupling  [16],  and  concentration  dependent  elastic 
modulus  2,17],  are  neglected  because  the  focus  here  is  to  investi¬ 
gate  the  impacts  of  plastic  yield  of  current  collector.  Therefore,  the 
diffusion  is  described  by 


9c  _  a 2_c 

at  “  a? 


(i) 


where  c  is  the  molar  concentration  of  Li-ions  in  the  plate,  D  is  the 
diffusivity  of  lithium  ions. 

Consider  a  full  cycle,  see  Fig.  3.  Both  charge  and  discharge  start 
with  galvanostatic  operation  followed  by  potentiostatic  operation. 
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Fig.  2.  The  symmetric  layered  electrode  investigated  in  this  work. 


In  the  charge  phase  of  the  cycle  the  lithiation  of  anode  starts  from  a 
lithium  free  state.  Hence,  the  initial  condition  and  the  boundary 
conditions  for  galvanostatic  charging  are 

c  =  0  for  t  =  0  (2a) 


d-  =  In¬ 
ez  F 


for  z  =  h\ 


(2b) 


D—  =  0  for  z  =  0  (2c) 

oz 

where  F  =  96485.3  C  mol-1  is  Faraday’s  constant  and  in  is  the 
surface  current  density  which  is  positive  for  lithiation  but  negative 
for  delithiation.  Due  to  the  symmetry,  only  the  equations  of  the 
upper  active  layer  are  provided.  The  distribution  of  Li-ion  concen¬ 
tration  is  provided  by  Crank  [18]: 
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It  is  assumed  that  the  charge  operation  switches  from  the  gal¬ 
vanostatic  to  the  potentiostatic  at  the  time  moment  to  with  Li-ion 
concentration  c(z,to).  Hence,  the  initial  and  boundary  conditions 
for  potentiostatic  charging  are 


c  =  c(z,  to)  for  t  =  t0 

(4a) 

c  =  c0  for  z  =  h i 

(4b) 

Fig.  1.  Explanation  of  the  terms  “pure  elastic  deformation”,  “plastic  shakedown”  and  “cyclic  plasticity”  with  the  stress-strain  curve  of  a  strain  controlled  uniaxial  tension/ 
compression  test.  aY  is  the  yield  strength. 
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Fig.  3.  Illustration  of  one  cycle  of  charge  and  discharge. 


D—  =  0  for  z  =  0  (4c) 

where  Co  is  the  constant  concentration  at  the  charging  surface.  The 
expression  for  c(z,to)  is  obtained  by  Equation  (3)  with  insertion 
t  =  to-  The  solution  can  be  obtained  by  using  the  method  of  sepa¬ 
ration  of  variables  and  is  written  as 


where  eo  is  the  in-plane  strain  of  the  plate  electrode,  Q  is  the  partial 
molar  volume,  E  =  E/(  1  -  v)  is  the  biaxial  modulus  and  v  is  the 
Poisson’s  ratio. 

The  current  collector  is  considered  as  an  elastoplastic  material 
with  linear  kinematic  hardening  behavior,  see  Fig.  4.  Its  elastic 
modulus  is  Ec,  plastic  modulus  Ep,  and  yield  strength  oY.  In  charge 
process  the  current  collector  is  stretched  due  to  the  expansion  of 
active  materials,  while  in  discharging  the  current  collector  is 
unloaded  due  to  the  shrinkage  of  active  materials. 

When  the  current  collector  is  elastic,  the  biaxial  stress  gc  de¬ 
velops  linearly  with  the  in-plane  strain  eq.  In  the  linear  hardening 
stage  after  the  current  collector  yields,  the  equivalent  stress 

o e  =  ^3SySy/2  changes  linearly  with  the  equivalent  plastic  strain 

ep  =  ^2^.^./3.  In  the  elastic  unloading  stage,  the  incremental 

stress  Aac  is  linear  to  the  incremental  strain  Aeo-  In  the  reversing 
hardening  stage  the  incremental  equivalent  stress  is  proportional 
to  the  incremental  equivalent  plastic  strain.  The  cyclic  plastic 
constitutive  behavior  is  depicted  in  Fig.  4  and  described  as 

oc  =  E'ce o  for  ep  =  0  (8a) 

A ae  =  EpAep  for  ep  ^  0,  dsp  >  0  (8b) 

Aac  =  E'Ae0  for  ep^ 0,  d£p  =  0  (8c) 
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In  the  discharging  phase  of  the  full  cycle  the  delithiation  will 
start  with  galvanostatic  operation  from  the  saturation  state  having 
homogenous  lithium-ion  distribution,  and  continue  with  poten- 
tiostatic  operation  from  the  time  moment  t  =  t2-  The  delithiated  Li- 
ions  of  galvanostatic  discharging  can  be  obtained  by  simply  using 
Equation  (3)  with  the  replacement  of  t  by  discharging  time  (t  -  ti). 
And  the  delithiated  Li-ions  during  the  follow-up  potentiostatic 
discharging  can  be  obtained  by  Equation  (5)  with  replacement  of  t 
and  to  by  discharging  time  (t  -  t\)  and  (t2  -  ti). 

Considering  that  the  layered  electrode  is  large  and  the  Li-ions 
are  inserted  symmetrically,  the  electrode  does  not  bend.  Among 
the  totally  six  stress  components,  only  the  two  in-plane  biaxial 
stresses  gx  and  ay  are  non-zero.  The  out-of-plane  stress  gz,  which 
arises  from  the  compression  by  adjacent  electrodes,  is  neglected  as 
it  is  small  and  considered  incapable  of  significantly  impacting  the 
elastoplastic  behavior  of  current  collector.  Because  ax  =  ay,  we  use 
g\  and  crc  to  denote  the  biaxial  stress  in  the  active  layer  and  current 
collector,  respectively.  According  to  the  equilibrium  condition,  o\ 
and  ac  must  satisfy 


/ii  o 

2  J  G^dz+  J  Gcdz  =  0  (6) 

o  -he 

Focusing  on  the  active  materials  like  graphite  and  LiFeP04  which 
do  not  show  significant  plastic  behaviors,  the  active  material  layers 
are  assumed  to  be  isotropic  and  always  elastic.  In  the  active  layer 
the  stress  g\  is 

=  E'(£0-Qc/3)  (7) 


where  E'c  =  Ec/(  1  -  v)  is  the  biaxial  modulus  of  the  current  col¬ 
lector.  Note  that  Sy  =  Gy  -  Gkk^ij/3  is  the  deviatoric  stress,  and  the 
plastic  strain  ef  =  -2ejj  =  -2f?y  due  to  the  assumption  of  plastic 
incompressibility,  we  have 


Oe  —  Gc 

(9a) 

CN 

II 

CN 

II 

% 

(9b) 

Substitute  Equations  (7)  and  (8a)  into  Equation  (6),  the  biaxial 
stresses  in  the  current  collector  and  active  layer  in  elastic  stage  are 

(10a) 

(10b) 

where  E\  =  1  +  l/2Ec/Eihc/hi 

is  a  dimensionless  variable  deter- 

mined  by  the  modulus 

ratio  and  thickness  ratio, 

ph\ 

Q  =  /  cdz  is  the  SOC  (state  of  charge)  and  cs  is  the 

Jo 

saturation  concentration  at  stoichiometric  limit. 


Fig.  4.  Illustration  of  the  plastic  constitutive  behavior  of  the  current  collector.  A^yde  is 
the  plastic  strain  range. 
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In  the  plastic  hardening  stage  of  charge 


2EpQcs  _ 
^2 


(11a) 


<Ti  =  E\Qc - 


he  &Y 
hi  2 T2  + 


3  T2  ^ 


(lib) 


where  r2  =  hc/h^Ep/E\  + 1  +  2EP/E'C.  In  the  elastic  unloading 
stage,  taking  the  end  of  charging  as  the  reference  state,  the  incre¬ 
mental  stresses  are 


Affc  =  Q  (12a) 

Ji  1 

Aff,  =  ^^AQ-l^GAc  (12b) 

In  this  stage,  reversing  yield  would  be  induced  if  Acrc  exceeds 
2 ay.  Whereas  if  Acrc  is  smaller  than  2oy,  there  would  be  no  reversing 
yield  and  the  current  collector  is  unloaded  and  reloaded  elastically 
in  the  following  cycles. 

In  the  reversing  hardening  stage  after  the  reversing  yield,  the 
incremental  equivalent  stress  A ae  is  again  linear  to  the  incremental 
equivalent  plastic  strain  be?.  Similar  to  Equations  (11a)  and  (lib), 
we  have 


A  (7  c 


2Epflcs 
3  A  KL 


(13a) 


A*!  =  (l+f)  ^AQ-^Ac  (13b) 

Finally  in  the  elastic  reloading  stage,  the  stresses  develop  again 
linearly  with  the  in-plane  strain  as  described  by  Equations  (12a) 
and  (12b). 

In  the  cycles,  the  plastic  strain  range  A^yde  is  very  important  as 
it  directly  impacts  the  low-cycle  fatigue  life  (Here  A^yde  is  used  to 
distinguish  from  the  incremental  plastic  strain  A^).  According  to 
Fig.  4,  the  plastic  strain  range  A^yde  is  equal  to  the  incremental 
plastic  strain  in  the  reversing  hardening  stage.  Therefore,  according 
to  Equations  (13a)  and  (8b) 

Kcle  =  ^^~2Qy)  (14) 

where  Qs  is  the  saturation  SOC  when  delithiation  starts,  Qy  is  the 
SOC  corresponding  to  the  initial  yield  and  can  be  determined  by 
Equation  (10a),  thus  Qs  -  2Qy  corresponds  to  the  change  of  SOC  in 
reversing  hardening  stage. 


3.  Results  and  discussions 


3.1.  Yield  criterion ,  plastic  shakedown  and  cyclic  plasticity 

Firstly  we  discuss  the  yield  criterion.  The  current  collector  yields 
when  the  equivalent  stress  reaches  yield  strength,  i.e.  oe  =  ay. 
Introducing  Equations  (9a)  and  (10a),  we  have  the  yield  criterion 

jTTxa  -  (,5> 

TT  +  2  hj 

where  gy  =  3 gy/E\Qcs  is  the  dimensionless  yield  stress.  This  cri¬ 
terion,  depicted  in  Fig.  5  using  the  red  line  (in  the  web  version), 
indicates  that  active  materials  with  lower  elastic  modulus  E\  and 


Fig.  5.  Yield  map  showing  the  separation  of  three  types  of  elastoplastic  behaviors  of 
current  collector,  i.e.  pure  elastic  deformation,  plastic  shakedown  and  cyclic  plasticity, 
with  respect  to  the  elastic  modulus  ratio  Ec/Ei  and  thickness  ratio  hc/hi-  This  map  is  for 
fully  charged  graphite  electrodes. 


smaller  lithiation  induced  strain  Qcs  are  favored  to  avoid  the  plastic 
yield  of  current  collector.  In  addition,  the  elastic  modulus  ratio  and 
thickness  ratio  of  current  collector  to  active  material  also  impact 
the  yield. 

Once  yields,  the  current  collector  may  exhibit  two  types  of 
plastic  behaviors  in  the  following  cycles,  i.e.  plastic  shakedown  and 
cyclic  plasticity.  In  the  former  case  the  unloading  and  reloading  is 
purely  elastic,  see  Fig.  1(b).  While  in  the  latter  case  alternating 
plastic  yield  takes  place  in  each  loading-unloading  cycle,  see 
Fig.  1(c).  Consider  a  discharge  operation  from  saturation  state  to 
lithium-free  state  in  which  AQ  =  Qs,  reversing  yield  takes  place 
once  the  incremental  stress  A gc  exceeds  2gy.  Introducing  this 
condition  into  Equation  (12a),  the  critical  condition  distinguishing 
between  plastic  shakedown  and  cyclic  plasticity  is  obtained  as 


Tc 


1  Qs 

,  1  hc  2 
+  2H7 


<7y 


(16) 


This  critical  condition  is  plotted  in  Fig.  5  using  the  blue  line. 
Comparing  Equations  (16)  with  Equation  (15)  and  noting  Qy  is  the 
SOC  corresponding  to  the  first  yield,  it  is  clear  that  plastic  shake- 
down  takes  place  if  the  initial  yield  of  current  collector  occurs  after 
the  electrode  has  been  charged  half  full,  i.e.  Qy  >  0.5Qs,  whereas 
cyclic  plasticity  will  be  induced  when  Qy  <  0.50s. 

A  yield  map  is  provided  based  on  the  above  critical  conditions, 
see  Fig.  5.  This  figure  shows  separation  of  the  three  possible 
elastoplastic  behaviors  of  current  collector  in  a  fully  charged 
graphite  electrode.  It  shows  that  the  elastoplastic  type  changes 
from  pure  elastic  deformation  to  plastic  shakedown  and  further 
to  cyclic  plasticity  with  increasing  elastic  modulus  ratio  Ec/Ei  as 
well  as  decreasing  thickness  ratio  hclh\.  The  reason  is  that  the 
current  collector  with  higher  Ec  applies  stronger  restriction  to  the 
active  layer,  and  smaller  thickness  hc  results  in  smaller  cross- 
sectional  area  of  the  current  collector,  both  of  which  lead  to 
higher  stress  and  larger  plastic  strain  in  the  current  collector 
during  the  cycles. 

According  to  Fig.  5,  the  thickness  ratio  hc/hi  is  preferably  higher 
than  2.4  to  avoid  the  plastic  yield  of  current  collector  if  we  choose  a 
copper-graphite  electrode  whose  modulus  ratio  Ec/Ei  is  about  10. 
This  thickness  ratio  is  higher  than  that  employed  in  real  electrodes 
[19].  It  may  be  an  error  induced  by  assigning  the  maximum  lith¬ 
iation  induced  volumetric  strain  Qcs  =  8%  of  pure  graphite  to  the 
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graphite  composite  whose  Qcs  should  be  smaller  due  to  the  exis¬ 
tence  of  pores  and  other  non-active  compositions  in  the  composite. 
Therefore,  the  structural  designs  provided  by  this  paper  are 
considered  conservative. 

3.2.  Diffusion  induced  stresses 

In  the  following  sections,  impacts  of  the  plastic  yield  of  current 
collector  on  DISes  will  be  discussed.  The  current  collector  is  chosen 
as  copper  and  the  active  material  as  graphite.  The  employed  ma¬ 
terial  constants  of  copper  are  yield  strength  oY  =  300  MPa,  elastic 
modulus  Ec  =  150  GPa  and  plastic  modulus  Ep  =  3  GPa.  The  elastic 
modulus  of  graphite  is  E\  =15  GPa.  The  operation  of  charge  and 
discharge  follows  that  provided  by  Fig.  3. 

Fig.  6(a)  and  (b)  shows  the  stresses  along  plate  thickness  in  a 
time  moment  after  the  current  collector  has  yielded  in  charge  and 
discharge,  respectively.  The  thickness  ratio  hc/hi  is  set  to  0.25  ac¬ 
cording  to  Fig.  5  to  ensure  the  plastic  yield  of  current  collector.  It  is 
found  that  the  estimations  of  stresses,  in  both  current  collector  and 
active  layer,  are  significantly  higher  once  the  plastic  yield  of  current 
collector  is  neglected.  The  reason  is  that  a  substantial  amount  of 
diffusion  induced  stress  in  the  layered  electrode  arises  from  the 
constraint  of  active  layer  by  current  collector.  Once  the  current 
collector  yields,  the  plastic  relaxation  enables  the  layered  electrode 
occur  relatively  large  deformation  with  respect  to  small  increment 
of  stress.  Therefore,  ignoring  the  plastic  yield  of  current  collector 


Fig.  6.  The  dimensionless  biaxial  stress  a  =  3 o/E\Qcs  in  the  electrode:  (a)  in  charge 
process,  and  (b)  in  discharge  process. 


may  lead  to  overestimations  of  the  stresses  and  underestimations 
of  the  strength  of  electrodes. 

As  seen,  the  diffusion  induced  stresses  in  active  layers  can  be 
reduced  by  the  plastic  relaxation  of  current  collector.  Therefore, 
Fig.  7  is  further  prepared  to  illustrate  the  impacts  of  the  plastic  yield 
by  plotting  the  evolutions  of  the  DISes  in  certain  positions  of  the 
electrode  in  cycles.  According  to  the  yield  map  provided  by  Fig.  5, 


0.0  0.2  0.4  0.6  0.8  1.0 

soc 


Fig.  7.  Evolutions  of  the  biaxial  stress  against  SOC  in  cycles  of  charge  and  discharge:  (a) 
in  the  current  collector,  (b)  at  the  interface  between  current  collector  and  active  layer 
where  z/hi  =  0,  and  (c)  at  the  inlet  surface  where  z/h-i  =  1.  Three  thickness  ratios  are 
employed  according  to  Fig.  4  so  that  the  blue,  red  and  dark  lines  describe  the  cases  of 
pure  elastic  deformation  (PED),  plastic  shakedown  (PS)  and  cyclic  plasticity  (CP), 
respectively.  The  thin  lines  represent  the  stresses  in  the  first  cycle.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 
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the  thickness  ratio  hc//ii  is  set  to  3, 1.5  and  0.25,  respectively,  cor¬ 
responding  to  the  case  of  pure  elastic  deformation,  plastic  shake- 
down  and  cyclic  plasticity. 

Fig.  7(a)  shows  the  evolutions  of  the  stress  in  the  current  col¬ 
lector.  It  is  clearly  seen  the  change  from  pure  elastic  deformation  to 
plastic  shakedown  and  finally  to  cyclic  plasticity  with  respect  to 
decreasing  thickness  ratio.  In  the  case  of  hc/hi  =  3,  the  current 
collector  is  thick  so  that  it  cannot  be  stretched  to  yield,  the  stress  in 
current  collector  changes  linearly  with  SOC.  In  the  case  of  shake- 
down  when  hc/hi  =  1.5,  the  plastic  strain  is  limited  and  the  stress  in 
current  collector  changes  linearly  with  SOC  since  the  second  cycle. 
Nonlinear  behavior  occurs  only  in  the  first  charge.  In  the  case  of 
cyclic  plasticity  when  hc/hi  =  0.25,  the  current  collector  is  too  thin 
so  that  the  current  collector  yields  in  every  charge  and  discharge 
and  the  stress  evolves  in  the  form  of  a  closed  loop.  With  the  change 
of  plastic  type,  the  stress  in  current  collector  evolves  in  larger 
ranges  with  higher  peak  values  and  higher  residual  stresses. 

Fig.  7(b)  and  (c)  illustrates  the  stresses  in  the  position  of  z/fti  =  0 
and  z\h\  =  1  in  the  active  layer,  respectively.  It  is  seen  that  both  the 
peak  compressive  stress  and  cyclic  stress  range  in  the  charge/ 
discharge  cycle  are  significantly  decreased  when  the  elastoplastic 
type  changes  from  pure  elastic  deformation  to  cyclic  plasticity, 
although  tensile  stress  is  increased  due  to  the  residual  deformation 
of  current  collector.  This  figure  for  one  hand  reveals  one  mecha¬ 
nism  how  tensile  stress  is  generated  in  anode,  for  the  other  hand  it 
is  interesting  to  see  that  applying  thinner  current  collectors  and 
allowing  them  to  yield  is  beneficial  to  the  integrity  of  active  layer 
as  the  peak  stress  and  the  cyclic  stress  range  are  significantly 
decreased. 

The  designs  allowing  plastic  yield  show  two  advantages 
compared  with  that  free  of  plastic  yield.  Firstly,  the  thickness  of 
current  collector  is  decreased  significantly,  with  which  electrode 
capacity  can  be  enhanced  as  more  space  can  be  provided  for  the 
active  materials.  It  is  found  in  Fig.  7  that  the  thickness  of  current 
collector  in  the  case  of  plastic  shakedown  is  only  about  50%  of  that 
in  the  case  of  pure  elastic  deformation.  This  ratio  is  even  lower  in 
the  case  of  cyclic  plasticity.  Secondly,  the  cycle  life  of  electrode  may 
benefit  from  the  plastic  yield  because  both  the  peak  stress  and  the 
cyclic  stress  range  in  active  layers,  which  is  critical  to  the  me¬ 
chanical  fade  and  electrochemical  stability,  are  significantly 
decreased.  Therefore,  applying  a  thinner  current  collector  and 
allowing  it  to  yield  is  beneficial  not  only  to  the  capacity  but  also  to 
the  structural  integrity  and  cycle  life  of  electrode. 

It  should  be  noted  that  the  stress  in  active  layer  is  decreased 
with  the  sacrifice  of  higher  stress  and  cyclic  stress  range  in  the  thin 
current  collector.  However,  this  increase  of  stress  may  be  not  so 
important  to  the  battery  performance  because  current  collector, 
which  functions  as  a  passageway  of  electrons,  does  not  involve  in 
the  electrochemical  reaction  and  storage  of  Li-ions.  Therefore,  this 
increase  of  stress  in  current  collector  may  be  acceptable  on  the 
condition  that  the  plastic  deformation  does  not  lead  to  crack. 

3.3.  On  the  plastic  yield 

In  this  section  we  discuss  whether  the  plastic  yield  of  current 
collector  is  allowable.  If  yes,  which  type  is  better  between  plastic 
shakedown  or  cyclic  plasticity?  How  to  determine  the  thickness  of 
current  collector  in  designing? 

For  the  first  question,  it  is  known  plastic  yield  should  be  avoided 
for  ordinary  load-bearing  structures  because  cyclic  fatigue  may  lead 
to  structural  failure.  However,  the  electrode  of  a  Li-ion  battery  may 
be  not  so  sensitive  to  fatigue  for  two  reasons.  Firstly,  the  cycle  life  of 
a  real  battery  may  be  less  than  one  thousand  of  times  and  one  cycle 
of  charge  and  discharge  may  take  several  hours  to  several  days. 
They  are  considered  not  adequate  to  induce  fracture  in  the  case  of 


finite  plastic  strain  because  experiments  showed  that  poly¬ 
crystalline  copper  took  more  than  105  cycles  to  fracture  with  con¬ 
stant  and  low  plastic  strain  amplitude  [20].  Secondly,  the  materials 
of  current  collector  like  copper  and  aluminum  are  known  with  good 
ductility.  The  reported  ultimate  fracture  strain  of  copper  film  is  6% 
[21],  whereas  the  extreme  in-plane  strain  e0  of  the  copper- 
graphite  electrode  is  calculated  as  about  2.6%,  which  is  limited  by 
the  expansion  of  graphite  and  far  below  the  reported  value. 
Therefore,  finite  plastic  deformation  of  current  collector  may  be 
allowable. 

The  second  question  is  which  type  is  better  between  plastic 
shakedown  and  cyclic  plasticity.  Generally,  the  case  of  plastic 
shakedown  shows  good  balance  among  all  aspects,  e.g.  capacity 
enhancing  due  to  thinner  current  collector,  better  electrochemical 
stability  due  to  lower  stress  in  active  layer,  and  good  fatigue 
resistance  due  to  the  zero  plastic  strain  range  since  the  second 
cycle.  Electrode  structure  can  be  designed  according  to  the  critical 
criterion  of  plastic  shakedown  provided  by  Equation  (16).  Accord¬ 
ing  to  calculation,  changing  the  thickness  ratio  hc/hi  from  2.4  to  1.1, 
corresponding  to  the  transition  from  pure  elastic  deformation  to 
plastic  shakedown,  may  lead  to  a  reduction  of  the  volume  of  cur¬ 
rent  collector  by  more  than  50%  as  well  as  the  decrease  of  stress  in 
active  layer  by  40%. 

Compared  with  plastic  shakedown,  the  case  of  cyclic  plasticity  is 
even  superior  in  the  aspects  of  capacity  enhancing  and  electro¬ 
chemical  stability,  though  the  plastic  yield  in  every  reversal  may 
lead  to  low-cycle  fatigue  and  significant  shortening  of  cycle  life.  The 
number  of  cycle  to  failure  in  low-cycle  fatigue  is  usually  charac¬ 
terized  by  the  Coffin-Manson  relation  [22,23]: 

Kde  =  ^  07) 

where  N  is  the  expected  cycle  number  to  failure,  e'j  is  the  fatigue 
ductility  coefficient  and  c  is  the  fatigue  ductility  exponent.  Being 
aware  that  both  and  c  are  known  empirical  constants  for  a  given 
material,  the  allowed  plastic  strain  range  A^yde  thus  can  be 
determined  according  to  the  designed  cycle  numbers  N.  Finally,  the 
thickness  ratio  hc/hi  can  be  looked  up  from  Fig.  8.  This  designing 
method  enhances  electrode  capacity  and  decreases  the  stress  in 
active  layer  as  much  as  possible  with  the  aim  of  an  expected  cycle 
life. 


Fig.  8.  The  in-plane  strain  e0  of  graphite-copper  electrode  and  the  plastic  strain  range 
A^cie  at  saturation  state  vs.  thickness  ratio  hclh\.  The  peak  value  of  in-plane  strain  is 
about  2.6%  when  hc  =  0,  consistent  with  the  lithiation  induced  volumetric  strain  8%  of 
free  expansion  graphite. 
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4.  Conclusions 

This  work  studies  the  impacts  of  plastic  yield  of  current  collector 
on  the  diffusion  induced  stress  in  a  symmetric  layered  electrode. 
Based  on  analytical  formulations,  three  types  of  elastoplastic  be¬ 
haviors  of  current  collector,  i.e.  pure  elastic  deformation,  plastic 
shakedown  and  cyclic  plasticity,  are  identified.  Criterions  sepa¬ 
rating  the  three  cases  are  proposed.  It  is  found  the  plastic  behavior 
of  current  collector  would  change  from  pure  elasticity  to  plastic 
shakedown  and  finally  to  cyclic  plasticity  with  increasing  elastic 
modulus  ratio  Ec/Ei  and  decreasing  thickness  ratio  hc/hi  of  current 
collector  to  active  layer. 

As  for  the  diffusion  induced  stress,  it  is  found  that  neglecting  the 
plastic  yield  of  current  collector  would  result  in  overestimations  of 
the  stress  in  whole  electrode  and  underestimation  of  electrode 
strength.  By  simulating  the  evolution  of  stress  in  cycles  of  charge 
and  discharge,  it  is  suggested  that  applying  a  thin  current  collector 
and  allowing  finite  plastic  deformation  is  beneficial  not  only  to  the 
battery  capacity  as  more  space  can  be  provided  for  active  materials, 
but  also  to  the  electrochemical  stability  and  cycle  life  of  electrode 
because  the  stress  in  active  layer  is  significantly  decreased. 

Regarding  the  plasticity  of  current  collector,  firstly  finite  plastic 
deformation  is  considered  acceptable  due  to  the  few  cycle  numbers, 
long  cycle  period  and  limited  plastic  strain  of  current  collector. 
Secondly,  electrode  structures  responsible  for  plastic  shakedown 
and  cyclic  plasticity  are  discussed  and  compared.  The  case  of  plastic 
shakedown  shows  good  balance  among  all  aspects,  e.g.  capacity 
enhancing  due  to  thinner  current  collector,  better  electrochemical 
stability  due  to  lower  stress  in  active  layer,  and  good  fatigue 
resistance  due  to  the  zero  plastic  strain  range.  In  contrast,  the  case 
of  cyclic  plasticity  is  even  superior  in  the  former  two  aspects  but 
with  the  sacrifice  that  the  cycle  life  is  significantly  shortened  due  to 
low-cycle  fatigue.  Therefore,  structural  design  is  provided  by  the 
critical  criterion  of  plastic  shakedown  in  the  former  case  while  by 
the  Coffin— Manson  relation  according  to  an  expected  cycle  life  in 
the  latter  case. 

It  should  be  noted  this  paper  discusses  the  impact  of  plastic 
yield  by  taking  a  graphite-copper  anode  for  example.  However,  the 
main  conclusion  of  this  work,  i.e.  applying  a  thin  current  collector 
and  allowing  it  to  yield  is  beneficial  to  both  the  capacity  and  cycle 
life,  is  also  applicable  for  cathode.  The  only  difference  is  the  sign  of 
stress.  For  a  cathode  which  is  initially  lithium-saturated  and 
deformation  free,  the  current  collector  would  be  compressed  to 
yield  in  charge  (delithiation)  operation.  In  addition,  applying  a  thin 


current  collector  in  cathode  and  allowing  it  to  yield  results  in 
decrease  of  tensile  stress  in  active  layer. 
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